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Abstract: We propose a set of lattice measurements which could test whether the decon- 
fined, quark-gluon plasma, phase of QCD shows strong coupling aspects at temperatures 
a few times the critical temperature for deconfinement, in the region where the conformal 
anomaly becomes unimportant. The measurements refer to twist-two operators which are 
not protected by symmetries and which in a strong-coupling scenario would develop large, 
negative, anomalous dimensions, resulting in a strong suppression of the respective lattice 
expectation values in the continuum limit. Special emphasis is put on the respective opera- 
tor with lowest spin (the spin-2 operator orthogonal to the energy-momentum tensor within 
the renormalization flow) and on the case of quenched QCD, where this operator is known 
for arbitrary values of the coupling: this is the quark energy-momentum tensor. The pro- 
posed lattice measurements could also test whether the plasma constituents are pointlike (as 
expected at weak coupling), or not. 
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1. Introduction 

The heavy ion experiments at RHIC have given two rather surprising and important results, 
namely the medium effects know as elliptic flow and jet quenching turned out to be much 
larger than simple expectations based on perturbative QCD. This has led to the picture that 
the deconfined matter produced at RHIC is a nearly perfect fluid, so like a strongly coupled 
plasma (see, e.g., the review papers ]l], § and references therein). The coupling constant 
a s = g 2 /4-7r in QCD can never become large, because of asymptotic freedom, but it can be 
of order one at scales of order Aqcd ; and this might lead to an effectively strong-coupling 
behaviour. It is notoriously difficult to do reliable estimates in QCD when a s ~ 1, so it 
has become common practice to look to the strongly coupled J\f = 4 supersymmetric Yang- 
Mills (SYM) theory for guidance as to general properties of strongly coupled field theories 
at finite temperature (see the review papers ||, || for details and more references). Since 
conformal symmetry is an essential property of M = 4 SYM, this theory is probably not a 
good model for the dynamics in QCD in the vicinity of the deconfinement phase transition, 
where the conformal anomaly associated with the running of the coupling in QCD is known 
to be important. But lattice studies || show that the relative conformal anomaly (e — 3p)/e 
(e is the energy density and p is the pressure) decreases very fast with increasing T above T c 
and becomes unimportant (smaller than 10%) for temperatures T > 2T C ~ 400 MeV. Hence, 
there is a hope that, within the intermediate range of temperatures at 2T C < T < 5T C , which 
is the relevant range for the phenomenology of heavy ion collisions at RHIC and LHC, the 
dynamics in QCD may be at least qualitatively understood by analogy with N = 4 SYM 
theory at strong coupling. 

One result suggested by SYM theory is that all the leading twist operators that occur 
in the operator product expansion of deep inelastic scattering, with the exception of the 
energy-momentum tensor, should have vanishing expectation values in a strongly coupled 
plasma. For the N = 4 SYM plasma, whose strong-coupling limit can be studied via the 
gauge/gravity duality @, |§ , this result follows from the fact that only protected operators 
- those whose anomalous dimension is zero because of a symmetry, or conservation law - 
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do not acquire large negative anomalous dimensions [0, ||, |9|, [To|, 11 , [T^] . This result is in fact 
natural in any field theory whose coupling is large. As one measures this theory at smaller 
and smaller space-time scales, one uncovers more and more strong evolution (branching) of 



the quanta of the theory O, 14, 15, 16]. The smallness of the higher dimensional operators 



in the leading-twist series just represents the fact that the higher energy-moments of "bare" 
quanta are naturally small at strong coupling, because the energy has been shared among 
many quanta via the branching process. At finite temperature, it is natural to assume that 
the branchings have taken place between the temperature scale T and the "hard" resolution 
scale Q, with Q 3> T, at which the operator is evaluated. 

It is important to emphasize that the renormalization flow of the operators and, in par- 
ticular, their anomalous dimensions are determined by the vacuum properties of the theory 
- the temperature enters only as the natural scale at which this flow should begin (and 
which therefore controls the early running of the coupling in a theory like QCD). In partic- 
ular, at weak coupling, the anomalous dimensions are computable in the zero-temperature 
perturbative expansion, which is a series in powers of a s . This should be contrasted with the 
calculation of thermal expectation values, so like the pressure, whose perturbative expansion 
is quite subtle already at weak coupling, and in particular non-analytic in a s , because of 
infrared problems associated with the thermal Bose-Einstein distribution p7| , |l8| , Thus 
a non-perturbative study of the renormalization flow of the (unprotected) leading twist oper- 
ators would allow one to distinguish between genuine strong-coupling effects and the failure 
of the perturbation theory due to medium effects, which occurs already at weak coupling. 
This would avoid the ambiguity inherent in the present lattice studies of the QCD thermo- 
dynamics §, whose results cannot be accommodated by a strict perturbative expansion, yet 
they appear to be consistent, for temperatures T > 2.5T C , with the predictions of the HTL- 
resummed perturbation theory |2C, 21, 22, p3|, and not too far away from the strong-coupling 
limit of M = 4 SYM H, [2§. 

In particular, recent lattice calculations [^] of the fluctuations of the electric charge, 
baryonic number, and strangeness in the quark-gluon plasma appear to be remarkably close 
to the respective results of HTL-resummed perturbation theory |27j (and also to the ideal 
gas limit) already for T > 1.5T C , thus strongly supporting a quasiparticle picture of the weak 
coupling type. The fact that the approach towards the ideal gas limit when increasing T above 
T c appears to be faster for the quark susceptibilities than for the pressure or energy density, 
is perhaps to be attributed to the fact that the conformal anomaly, which is so important for 
the thermodynamics in the vicinity of T c , is less important for fermionic observables, like the 
above susceptibilities. 

Now, there is a priori no contradiction in having a quasiparticle picture also at strong 
coupling, as shown by the fact that the entropy density of the N = 4 SYM plasma in the 
strong-coupling limit is close to the respective value at zero coupling. However if the effective 
coupling is large, one expects the quasiparticles to be highly composite, without a pointlike 
core carrying a significant fraction of the quasiparticle energy. This is illustrated by recent 
calculations of deep inelastic scattering (DIS) off the strongly coupled N = 4 SYM plasma, 
which show that, when this plasma is measured on a hard resolution scale Q ^> T, one finds 
only low-energy constituents, with energy fractions x < T/Q < 1 [15, [l6|]. The higher Q is, 
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the smaller are the energy fractions, meaning that there are no pointlike constituents. 

How to determine what is the corresponding picture for the quark-gluon plasma ? Of 
course, one cannot literally perform a deep inelastic scattering on the QCD matter produced 
at RHIC to test whether or not there are pointlike constituents having energy of order of 
T (in the rest frame of the plasma). The only experimental evidence on this comes from 
the phenomenology of "jet quenching", which within perturbative QCD at least, is the pro- 
cess responsible for both energy loss and transverse momentum broadening of a hard probe 
propagating through the medium. The relevant transport coefficient q (the "jet-quenching 
parameter") is given by [p8| 

. _ 4TT 2 a s N c dxG(x,Q 2 ) 

q N 2 -l dV ' 1 j 

where dxG(x, Q 2 )/dV is the number of gluons per unit volume in the plasma measured on 
the relevant energy (x) and virtuality (Q 2 ) resolution scales. It is generally assumed that Q is 
of the order of the saturation momentum of the plasma, since this is the typical momentum of 
the gluons exchanged between the jet and the medium. Weak-coupling estimates of q using 
ideal gas formulas for the density of the plasma constituents at the scale T together with 
perturbative evolution to the hard scale Q give q ~ (0.5 -j- l)GeV 2 /fm, while phenomenology 



1 29, 30 ] rather suggests that q should be somehow larger, between 5 and 15 GeV /fm. This 



difference supports the picture of strong evolution in the plasma, and hence of strong coupling 



1 31 ] . One should nevertheless keep in mind that this phenomenology is quite difficult and not 
devoid of ambiguities: strong assumptions are necessary in order to compute q, and also to 
extract its value from the RHIC data (see, e.g., the discussion in |j32^ ). 

In view of the experimental difficulties, it is natural to ask whether lattice gauge the- 
ory can illuminate this question. Computing the DIS structure functions on the lattice is 
in principle possible: via the operator product expansion and for sufficiently high-Q 2 , the 
moments of the structure functions can be related to expectation values of operators with 
spin n and (classical) dimension n + 2 — the leading twist operators — which form an infinite 
series (only the even values of n being relevant for DIS). Given the space-like kinematics of 
the DIS process, these operator expectation values are effectively Euclidean, and thus can be 
evaluated on the lattice. In order to reconstruct the structure functions from their moments, 
one would need to measure a large number (in principle, infinite) of the latter, which is prac- 
tically tedious, if not impossible. Indeed, operators with spin n = 4, 6, 8, ... involve too many 
derivatives to be accurately evaluated in lattice QCD, although some attempts were done in 
that sense, for the case of the proton structure functions (see, e.g., [j33|). 

However, in order to answer the limited questions that we address here, such a full 
reconstruction of the plasma structure functions is actually not needed. What we instead 
propose is to measure the expectation value of the unique leading-twist operator with n = 2 
which is not protected by symmetries, and thus check whether the corresponding result is 
rapidly vanishing when approaching the continuum limit — as expected for a strong-coupling 
dynamics — , or rather it is slowly evolving away from the respective ideal gas expectations 
— as it should be the case at weak coupling. 
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Specifically, consider the two leading-twist operators with n = 2 in QCD, that is 

Of = q^iD v q - (trace) , (1.2) 

and, respectively, 

Of = -FIT FT + . (1.3) 

(It is understood that the fermionic operator Of involves a sum over quark flavors and a 
symmetrization of the Lorentz indices, and we neglect the masses of the quarks.) These 
two operators are well defined only with a renormalization prescription, and thus implicitly 
depend upon the resolution scale Q 2 . Since they have the same quantum numbers, they mix 
with each other under the renormalization flow. The following linear combination yields the 
total energy-momentum tensor, 

which is a conserved quantity, and thus is insensitive to quantum evolution (it does not depend 
upon Q 2 ). Clearly, this operator cannot be used to test whether the plasma has pointlike 
constituents, or not. Within perturbation theory, it is always possible to construct the linear 
combination of Of and Of which is orthogonal to T^ v within the renormalization flow 
and therefore vanishes in the continuum limit Q 2 — > oo (the respective anomalous dimension 
being negative). This is the operator whose expectation value we would like to measure on 
the lattice. But if the coupling is strong, we do not know how to explicitly construct this 
orthogonal combination. 

Fortunately, there is a simpler version of the theory where the identification of this 
operator becomes possible for any value of the coupling: this is quenched QCD. Loosely 
speaking, this is the theory obtained from QCD after removing all the quark loops. On the 
lattice, this is non-perturbatively defined by removing the fermionic determinant from the 
QCD action. Note that the quark fields are still present in this theory, but only as external 
probes. In particular, it makes sense to evaluate the fermionic operator (|l.2| ) in quenched 
QCD: at finite temperature, this amounts to computing the Matsubara Dirac propagator in 
the background of the thermal fluctuations of the gauge fields. Such a calculation effectively 
resums all the respective Feynman graphs of QCD, except for those involving quark loops. 
For instance, if the coupling is weak (i.e., for high enough temperatures) and for a given 
resolution scale Q 2 which is not too hard, the expectation value {Of (Q 2 ))t should be close 
to the respective value for an ideal fermionic gas, as given by the Fermi-Dirac distribution. 
Moreover, when Q ~ T, the temperature T is the only scale in the problem, so by dimensional 
arguments we expect 

(Of(Q 2 ~ T 2 ))t oc T 4 for any value of the coupling. (1-5) 

(We implicitly assume here that the temperature is sufficiently high for the QCD trace 
anomaly to become unimportant; in practice, T > 2T C .) On the other hand, in the con- 
tinuum limit Q 2 — ► oo at fixed T, the above expectation value must vanish: 

(Of(Q 2 -> oo)) T -» (fixed T) . (1.6) 
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Eq. (|1.6| ) will be derived in Sect. 4, but it is easy to see how it comes about. The quark 
can emit gluons — the more so, the harder the scale at which one probes its substructure. 
But the emitted gluons, as well as those from the thermal bath, are not allowed to emit 
quark-antiquark pairs. Hence, when the system is probed on a sufficiently hard scale, most 
of the total energy appears in the gluon fields. We thus see that, within quenched QCD, Oj U 
is the n = 2 operator orthogonal to the (total) energy-momentum tensor. In the continuum 
limit, the latter reduces to its gluonic component: T^ v — ► Og V (Q 2 ) as Q 2 — ► oo. 

Whereas Eq. (|1.G| ) holds for any value of the coupling, the rapidity of the evolution 
with increasing Q 2 — i.e., the rate at which (Oj U )t approaches to zero — depends upon 
the strength of the interactions. For a weak coupling, this evolution would be quite slow; 
using lowest-order perturbative QCD, we shall estimate in Sect. 4 that for a temperature 
T ~ 3T C and an inverse lattice spacing a^ 1 = Q ~ 4 GeV, the deviation of (O® )t from the 
corresponding ideal gas value should not exceed 30%. On the other hand, if the evolution 
is more like at strong coupling, and if the measurements of q are indicative of what should 
be expected in the strongly-coupled QCD plasma, one would expect (O® )t to be reduced 
by a factor of 5 or more. Of course, all the conclusions that could be drawn in this way 
would strictly apply to quenched QCD alone. However, we expect real (unquenched) QCD 
to behave similarly (within the same range of temperatures), because the asymptotic freedom 
property of QCD is driven by gluon dynamics. 

2. Leading— twist operators: from weak to strong coupling 

Although the main emphasis in this paper is not on the process of deep inelastic scattering by 
itself, but rather on the lattice evaluation of specific, low spin, leading-twist operators, it is 
nevertheless natural to introduce these operators in the context of DIS and thus summarize 
some of their properties to be used later on. 

Within QCD, there are two infinite sequences of leading-twist operators: the fermionic 
ones, 

£>J n)w ""' Mn = g 7 {«(iD" 2 ) • • • (iD^)q - (traces) , (2.1) 
(the curly brackets around Minkowski indices mean symmetrization) , and the gluonic ones, 

O^ 1 '* = — - F^ lV (iD^ 2 ) • • • (iD^"- 1 )F^ n \ - (traces), (2.2) 
y 2 

where a trace over color indices is implicit. Such operators have spin n, classical dimension 
d = n + 2, and hence twist t = d — n = 2. For n = 2, we recover the operators expressing the 
energy-momentum tensor for quarks and gluons, respectively, cf. Eqs. (|1.2j)-( |0| ). For even 
values of n, n = 2, 4, 6, . . . , these operators enter the OPE of the current-current correlator 
which determines the cross-section for the standard DIS process, as mediated by the exchange 
of a space-like photon. More precisely, if the expectation values of the operators are evaluated 
directly at the resolution scale Q 2 for DIS (the virtuality of the space-like photon), then the 
OPE involves only the quark operators ( |2.1[) . In practice, however, it is convenient to evaluate 
the operators at some fixed renormalization scale fj,, or at some intrinsic physical scale - 
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say, the temperature for the case of DIS off the quark-gluon plasma. In such a case, the 
operators at the DIS scale Q 2 are obtained by following the renormalization flow from the 
original scale [i 2 (or T 2 ), and under this flow, the fermionic and gluonic operators having the 
same spin mix with each other. 

It is generally stated that the leading-twist operators dominate the OPE for DIS for 
sufficiently high Q 2 . This is strictly true only so long as the coupling is not too strong, as 
shown by the example of N = 4 SYM theory, where explicit calculations were also possible at 
strong coupling. To illustrate this, consider the leading-twist contributions to the moments 
of the DIS structure function F 2 (x, Q 2 ), which can be expressed as (we follow the conventions 
in Ref. §4|) 

1 dx x n ~ 2 F 2 (x,Q 2 ) ~ Af (Q 2 ) , (2.3) 

o 



where the approximate equality sign means that in the r.h.s. we have kept only the twist-2 

(n 

/ 



contribution. The quantity A^\q 2 ) is the expectation value of the spin-n fermionic operator 



(2.1) evaluated at the resolution scale Q 2 and with all the kinematical factors (responsible 
for the Minkowski tensor structure and for the actual dimension of the operator) stripped 
off. For instance, if F 2 refers to a proton with 4-momentum P^, then 

(P\0 ( f ^ n {Q 2 )\P) = A { f{Q 2 )2P^---P^ - (traces), (2.4) 

whereas for a plasma at temperature T : 

{0f ] ^^ n {Q 2 )) T = Af ] {Q 2 )T n 2n^ ■■■n^ - (traces), (2.5) 

where is the four-velocity of the plasma, with n M = (1,0,0,0) in the rest frame of the 
plasma. Note that A^ and F 2 are dimensionless in the case of the proton, but they have 
mass dimension two in the case of the plasma; this difference is related to the normalization 
of the proton wavefunction. The Bjorken x variable is defined as x = Q 2 /2(P • q) for DIS 
off the proton and, respectively, x = Q 2 /2T(n ■ q) for DIS off the plasma; here, q^ is the 
4-momentum of the virtual photon, with q^q^ = —Q 2 . 

Let us assume that the operator 0* is normalized at the scale an d ignore the issue of 
operator mixing for the time being (we shall return to this issue in Sect. 3). The corresponding 
operator at a different resolution scale Q is obtained by solving the renormalization group 
equation (henceforth the Minkowski indices will be kept implicit) 



=^ ) 0f ) => 0f\Q 2 )=^j^ 1 f\^0f\ti), (2.6) 

» 



where 7^ is the corresponding anomalous dimension, which in QCD depends upon the scale 
because of the running of the coupling. Clearly, a similar evolution equation applies for the 
expectation value A^\q 2 ) of the above operator. It turns out that the anomalous dimensions 
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are always negative, with the notable exception of the energy-momentum tensor ( |1.4[) , for 
which 7 = 0. Hence, A^(Q 2 ) — ► as Q 2 — > oo for any n > 4, whereas for n = 2 we have 

l 

dx ^(x, Q 2 ) — > const, as Q 2 — > oo , (2-7) 



o 

which is simply the statement of energy-momentum conservation. These results imply that, 
when increasing Q 2 , F2(x,Q 2 ) is increasing at small x, but decreasing at large x: the evolu- 
tion acts to decrease the average value of the energy fraction of the partons in the wavefunc- 
tion. This should be expected given the physical picture of the evolution in terms of parton 
branching, as described in the Introduction. The rate of the evolution towards zero for the 
unprotected operators, and also the weight of the small-x partons in the sum-rule Q2,7| ), are 
however quite different at weak and respectively strong coupling, as we now explain. 

Consider weak coupling first. To lowest order in perturbative QCD, the anomalous 
dimensions are obtained as (for a generic leading-twist operator O) 

lO\M> = ~ao — ; = - 7 1 7 2/A2 V ' 2 - 8 ) 

where ao is a positive number and in writing the second equality we have used the one-loop 
expression for the QCD running coupling, with 6q = (H-^c — 2A r j)/3. Then Eq. fl2.6|) implies 



of\Q 2 ) 



^(^q/Aqcd) 
ln(QVA| CD ) 



af/b 



Ofiti), (2-9) 



which shows that the approach towards zero with increasing Q 2 is merely logarithmic. Still 
at weak coupling, consider the case of a conformal field theory, so like M = A SYM, where 
the coupling a = g 2 /An is fixed; then, 70 = — aoia/An) is a fixed number of 0(a), and 



( f\Q 2 ) 



Q 2 



2 I a 



(n) _a_ 

0f\nl), (2.10) 



so that the evolution is typically faster than in QCD, since it is not slowed down by the 
decrease of the coupling with increasing Q 2 . But for both QCD and N = 4 SYM, the anoma- 
lous dimensions are small ~ 0(g 2 ) at weak coupling, so the leading-twist operators dominate 
indeed the moments of the DIS structure functions at high Q 2 : the corresponding contribu- 
tions from higher-twist operators are suppressed by inverse powers of Q 2 with exponents of 
0(1). 

Consider also the energy-momentum sum-rule (2.7) : although F2(x,Q 2 ) does rise at 



small x, as expected, pQCD predicts that this rise is rather mild, so that the integral in 
Eq. ( |2.7| ) is dominated by rather large values of x, of order 0.1. This is confirmed by the 
experimental data at HERA, which can be parameterized by a law F2(x,Q 2 ) ~ 1/x^® ) 
where the effective exponent X(Q 2 ) rises slowly with Q 2 , but it remains relatively small: 
\{Q 2 ) ~ 0.15 -=- 0.3. This expresses the fact that, at weak coupling, the branching proceeds 
via bremsstrahlung and favors the emission of small-x gluons, whose number grows very fast, 



-7- 



but which carry only a tiny fraction of the energy of their parent partons. Accordingly, most 
of the total energy remains in the "valence" degrees of freedom at large x. Since this is true 
for arbitrarily high Q 2 , it is clear that these valence constituents can be viewed as pointlike. 

What is the corresponding situation at strong coupling ? Since the respective results are 
not known for QCD, we focus on the N = 4 SYM theory, whose strong-coupling limit has been 
addressed via the gauge/string duality. By "strong coupling", we more precisely mean here 
the limit in which the gauge coupling is weak, g 2 <l, but the number of colors is sufficiently 
large, N c 3> 1, for the 't Hooft coupling to be large: A = g 2 N c » 1. (Recall that the 't Hooft 
coupling is the relevant coupling for perturbation theory at large N c .) Via the AdS/CFT 
correspondence, the (gluonic) leading-twist operators are mapped onto excited string states 
- closed strings which rotate in the AdS§ space-time geometry. By computing the energy 
spectrum for such states, one can deduce the quantum dimensions A(n) = n + 2 — of 
the dual operators 0^ n \ and thus extract their anomalous dimensions j( n \ One has thus 
found @ 0, |TJ 



and, respectively, 



7 (n) ~ A 1/4 for 1 < n < V\, (2.11) 



7 (n) ^ i n Jl for n > y 7 !. (2.12) 

27r vA 



That is, the anomalous dimensions are again negative (except, of course, for the protected 
energy-momentum tensor), and moreover they are extremely large: of C(A 1//4 ) for the opera- 
tors with lower spin. Via Eq. (|2.6|) , this implies that all the leading-twist operators with the 
exception of T^ v are strongly suppressed at high Q 2 , and hence they become irrelevant for 
DIS: the respective structure functions are rather controlled by T^ u together with protected 
higher-twist operators which have zero anomalous dimensions. 

The fact that the anomalous dimensions are so large at strong coupling means that the 
branching process is very fast and, as a result of it, all partons have fallen at small values 
of x. This is further confirmed by the fact that the anomalous dimensions ( p.ll[ )-( ^.12| ) rise 
with n, showing (via the moments ( |2.3| )) that the support of the structure function is now 
concentrated at small values of x. 

As an illustration of the situation at strong coupling, let us recall the results for DIS off 
the M = 4 SYM plasma in the strong-coupling limit A — > oo (or N c —* oo). In that limit, the 
anomalous dimensions for the non-protected leading-twist operators become infinite, while 
the higher-twist protected operators cannot contribute to the DIS cross-sections because of 



the energy-momentum conservation. Accordingly, the explicit calculation in Ref. [15] finds 
that there is no power-like tail in F2(x, Q 2 ) at high Q 2 . More interestingly, it also finds that 
there is an exponential tail, 

F 2 (x,Q 2 ) ~ xN 2 Q 2 exp{-c(Q/Q s (x))} for Q > Q s (x) = | (2.13) 

(c is a number), which reflects a tunneling process, reminiscent of the Schwinger mechanism: 
the highly-virtual (Q 3> Q s (x)) space-like current can decay into charged partons via a 
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tunnel effect induced by a uniform force ~ T 2 , which represents the action of the plasma on 
the dipole fluctuations of the current in this large-Ac limit. 

The exponential in Eq. ( |2.13| ) can be alternatively rewritten as exp{— c(x/x s (Q))} with 
x s(Q) = T/Q, showing that, for fixed Q T, the DIS structure function is essentially zero for 
any x larger than x s (Q) <C 1. This reflects the fact that, via successive branchings, all partons 
have fallen at small values of x. And, indeed, for sufficiently small values x < x s (Q) (or, 
equivalently, for low enough virtualities Q < Q s {x) at a given x), the exponential suppression 



goes away, and one finds [15] 



rp N 2/3 

F 2 (x,Q 2 ) ~ xN 2 c Q 2 (— J for Q < Q s (x) . (2.14) 



These estimates are such that the energy-momentum sum-rule (2.7) is saturated by the 
partons along the saturation line, i.e., those having x ~ x s (Q) : 

T 2 [ dxF 2 (x,Q 2 ) ~ T 2 x s F 2 {x s ,Q 2 ) ~ A 2 T 4 . (2.15) 
J o 

As also emphasized above, this sum-rule reproduces the right order of magnitude for the 
energy density of the strongly-coupled plasma: (T 00 )t ~ N 2 T 4 . One can similarly check 
that the higher moments with n > 4 are power suppressed at high Q 2 : 

j\xx n - 2 F 2 (x,Q 2 ) ~ <A C 2 Q 2 ~ N 2 Q 2 f£) . (2.16) 

3. Evolution of n = 2 operators in QCD for a generic coupling 

In this section we describe the evolution of the n = 2 leading-twist operators. We focus on 
the respective flavor-singlet operators, of which there are two: the quark (O'j ) and gluon 
(Og U ) energy-momentum tensors displayed in Eqs. ( |1.2D -( |L~3|) . Our emphasis will be on the 
mixing between these two operators under quantum evolution, leading to two orthogonal 
eigen-operators: one which is a priori known for any value of the coupling, since this is 
protected by energy-momentum conservation and hence it is scale-independent — this is, of 
course, the total energy-momentum tensor, T^ u = Oj U + Og V — , and the other one which 
is not protected and hence it depends upon the renormalization scale Q 2 . The latter, that 
we shall denote as Q tiV {Q 2 ), is explicitly known in QCD only for sufficiently high Q 2 , where 
perturbation theory can be used to compute the matrix of anomalous dimensions (see e.g. 
Ch. 18 in |]34|]). Here, we are rather interested in the situation at generic, and relatively 
strong, coupling, so our subsequent developments will be necessarily formal and incomplete: 
we shall try and use physical constraints and guidance from M = 4 SYM theory in such a 
way to characterize the mixing matrix and the structure of as well as we can without 
performing explicit calculations in QCD. 

In full generality, the relevant renormalization group equations can be written in matrix 
form as (from now on we shall omit the Lorentz indices) 
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which features the 2x2 anomalous dimension matrix 7(/i 2 ) of the n = 2 leading-twist 
operators. To lowest order in perturbation theory, the scale dependence of 7, as encoded 
in the running coupling a s (fi 2 ), factorizes out from the matrix structure. In that case it 
is convenient to pursue the analysis by diagonalizing the 7 matrix, since the corresponding 
eigenvectors are scale-independent (see e.g. Ref. p3|] ). However, such a simplification does 
not occur for the general case at hand. It is then preferable to consider the formal solution 
to Eq. ( |3.1D , as obtained by integrating this equation from the conventional renormalization 
scale (1q to the physically interesting scale Q 2 . The solution reads 

I o g (Q 2 ) \ _ ( M 99 M gf \ ( OM) \ ( o 2) 

\o f (Q 2 )) \M h M ff ) \OMj' 

where the evolution matrix M = (Mij), with i, j = g or /, can be compactly written as 

Q 2 

<1m2 7(^M, (3.3) 



, Or 

M(Q 2 ,fj, 2 ) = Pexpi J 

Mo 



M 2 



where the symbol P in the r.h.s. indicates the /^-ordering of the product of matrices in the 
series expansion of the exponential. 

We shall now argue that, still in full generality, the matrix M has only two independent 
components. This follows from energy-momentum conservation: the requirement that T = 
Of + Og be scale-independent, that is, 

O g (Q 2 ) + O f (Q 2 ) = O g (n 2 ) + Of^ 2 ) , (3.4) 

implies two constraints on the components of the matrix M (since the operators at the original 
scale ixq should be viewed as two independent quantities), leading to 

Mf 9 = l-Mgg, Mgf = 1~ Mff, (3.5) 

and therefore 1 



M = 99 }t . (3.6) 




So far we have not taken into account the fact that the coupling could be strong. Recall 
that our objective is to give a test of the idea that QCD is strongly coupled at a scale which is 
a few times the critical temperature for deconfinement T c . So, let us assume that the coupling 
is strong at the scale /io at which one starts the evolution. (This is the scale to be identified 
with the temperature T when the evolution takes place in the quark-gluon plasma phase.) 
Then we expect M to have an eigenvalue which is extremely small, nearly zero, corresponding 
to the fact that the "unprotected" operator B has a large and negative anomalous dimension, 



1 There are of course similar constraints on the anomalous dimension matrix, which there imply Y/ s = — j gg 
and j g f = —7//. This means that, at any scale (J?, the matrix has the left eigenvector (1,1) with 

eigenvalue 7t = 0. But the other eigenvector, orthogonal to T, is generally scale-dependent. 
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which is exponentiated by the evolution. Let us give a more formal argument in that sense: 
to that aim, we divide the logarithmic phase-space for the evolution \ii(Q 2 /^q) into a large 
number N of small steps with width e = (I / N) ln(Q 2 / fify , in such a way as to ensure that, 
within each interval, the anomalous dimension matrix is essentially constant. Then we can 
break the /i 2 '-ordered exponential in Eq. ( |3,3j ) into a product of N ordinary exponentials: 

M(Q 2 ,nl) = e^e 6 ™" 1 • • • e e71 , (3.7) 

where, of course, the quantities 7$ = 7(//f) are 2x2 matrices. The determinant det M is 
equal to the product of the determinants of the N matrices in the r.h.s. For any such a 
matrix, we can diagonalize 7, locally at /if : 7$ = hi diag(7r, 7e(i)) h^ 1 , where jt = (this is 
the anomalous dimension of the protected operator T), whereas 7e(i) is strictly negative (this 
is the anomalous dimension of the unprotected operator 0(/i 2 ) at fj, 2 = fj, 2 ). Then, clearly 



det e e71 = e e7e « det M = exp { / -^7e(// 2 ) > , (3 




where the integrand in the exponent is negative at any n 2 . Now, let us assume that for /i 2 
close to the lower limit /i 2 , (or anywhere else along the way from /j,q to Q 2 ), the anomalous 
dimension is extremely large, so like at strong coupling: this implies that det M ~ 0, as 
anticipated. By imposing detM = in Eq. ( |3.6| ), one finds Mff = 1 — M gg . We thus finally 
deduce the following, particularly simple, expression for the evolution matrix (with m = M gg ) 

M= f™ ™ V (3.9) 
\ 1 — m 1 — m I 

valid when the evolution takes place at least partially in a region in /j, 2 where the coupling is 
strong. Using this form for M in Eq. ([T^), one finds 

O g {Q 2 ) = mT, O f (Q 2 ) = (l-m)T, (3.10) 

which allows us to identify the n = 2 operator orthogonal to T, i.e., the one which has evolved 
essentially down to zero on the resolution scale Q 2 : 

G(Q 2 ) = {l-m)O g {Q 2 ) - mO f (Q 2 ) = 0. (3.11) 



If the quantity m = m(Q 2 , /Xq) were known theoretically, then Eq. ( 3.11 ) would be a prediction 
that could be tested in lattice gauge theory. Unfortunately, we do not know how to determine 
this quantity within the scenario that QCD is strongly coupled at low scales. 

What we do know, however, is that m should be independent of the precise value of the 
scale /Uq at which one starts the evolution: indeed, m is rather determined by the largest value 
of fi 2 at which the anomalous dimension 7©(/i 2 ) is still large. To see this, let us introduce an 
intermediate scale /i 2 , with < fi 2 < Q 2 , and thus write M = M-yM^, with 

1 

M\ = P exp ^ / T^-7(M 2 U al,d - ]> ^\'{ I ^'M'-)}- (3.12) 
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Now, assume that [i\ is such that the coupling is still strong in its neighborhood, so that M\ 
has the structure shown in Eq. ( |3,9|) with m — > mi. Then one can easily check that M = Mi, 
that is m = mi, and this even for a matrix M2 which has the most general possible structure, 
as shown in Eq. (|3,6|). (But of course in QCD we would also expect M2 to be of the simpler 
form ( |3.9| ) , since if the coupling is strong at some scale fi\ , it is still strong at the softer scale 
Ho < fJL±.) Hence, if fi$ is the largest value at which the coupling is still effectively strong, 
then we have m((5 2 ,/Ug) = m{Q 2 ,n 2 s ) for any [1q < us- 

Finally, one may worry that in QCD anomalous dimensions are scheme dependent and 
that there is no meaning to say that 7 is large. However, when Q 2 is large, the operators on the 
left hand side of Eq. fl3.2| ) have very little scheme dependence because a s {Q 2 ) becomes small 
at large Q 2 . The scheme dependence refers merely to the ability to transfer contributions 
between the evolution matrix M(Q 2 , [1q) and the operators O g {fjL 2 ) ) and Of{^) at the original 
scale. If QCD behaves like a strongly coupled field theory, then the operators O g (Q 2 ) and 
Of(Q 2 ) at the final scale are expressible in terms of the (protected) energy-momentum tensor, 
as shown in Eq. ( |3.10| ) . We have modeled our discussion to ressemble the situation in J\f = 4 
SYM theory (where there is no scheme dependence, because of the conformal symmetry), 
but we recognize that in QCD one could choose schemes in which a condition like Eq. ( |3.11| ) 
- i.e., the vanishing of O at the scale Q 2 — does not follow from the evolution (i.e., from 



the particular structure (3J3) of the evolution matrix M), but rather from the fact that a 
relation between O g and Of similar to ( 3.11 ) holds already at the original scale Hq. 



4. Evolution of n = 2 operators in quenched QCD 

Because we are unable to specify a definite value for the quantity m in Eqs. Q and flOU) , 
it is difficult to devise a test of strong coupling behaviour in terms of n = 2 leading-twist 
operators using lattice gauge theory for full (unquenched) QCD. However, experience with 
lattice calculations shows that there is generally not a large difference between quenched 
and unquenched QCD. Thus if full QCD is effectively a strongly coupled theory in the soft 
momentum region, one would naturally expect the same to be true for quenched QCD. As 
mentioned in the Introduction, quenched QCD consists in ignoring the quark loops, so the 
matrix element 7/ g of the anomalous dimension matrix in Eq. ( |3.1| ) must vanish. (Recall 
that this element describes a transition from gluon to quark fields.) Since j gg = —7/3 by 
energy-momentum conservation, and similarly jff = — 7 ff /, we deduce that the 7 matrix has 
a very simple structure in quenched QCD: 




t(^) = r ,JJ (4.1) 

V° vf J 

This structure is already telling us that the operator orthogonal to the energy-momentum 



tensor T = Of + O g is simply the quark operator Of. (Indeed, the matrix (4.1) has the 
left eigenvector (0,1) with eigenvalue jff < 0.) It is furthermore clear that the rest of the 
discussion of the renormalization group evolution for n = 2 goes exactly like in the previous 
section, so in particular Eqs. ( |3.S| ) and ( |3,11 ) are still true, but now with m = 1 (since gluon 
fields cannot change into fermions). Once again, Eq. ( [3.11 ) with m = 1 confirms that = Of. 
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Thus, within quenched QCD, a strong-coupling scenario predict Of(Q 2 ) ~ for a sufficiently 
hard scale Q 2 . At finite temperature, this in turn implies that the average value of the energy 
carried by a bare quark (one which is measured on a hard resolution scale Q 2 3> T 2 ) which 
is in equilibrium with a strongly-coupled thermal bath of gluons is very small, 



{q l0 iD q (Q 2 )) T * 0. 



(4.2) 



and in particular much smaller than the corresponding ideal-gas value (the Stefan-Boltzmann 
law for a gas of free, massless, quarks): 



(q 7 0iD q)P = NjNj^T* 



(4.3) 



By contrast, in a weak coupling scenario, the corresponding lattice result should be 
rather closed to the above ideal gas value, and slowly departure from it with decreasing 
lattice spacing a = 1/Q. One can easily evaluate the leading order perturbative corrections 
to fl4.3j ), and thus get a better estimate for what should be the result at weak coupling: using 



iff = — a ff 

in) 



a s (n 2 



a ff 



one finds (cf. Eq. ( |2.9[) with — >• 

(o f (Q 2 )) 
(oM)) 



4tt 13 3 

ajf and 6q = HiV c /3) 



Cf, 



(4.4) 



I^o/Aqcd) 
ln(Q 2 /A QCD ) 



8C F /36o 



(4.5) 



For example, for Q = 4 GeV, A QCD = 200 MeV, and fi = 3T C ~ 600 MeV, one finds that 
the perturbative evolution reduces the ideal-gas result (4.3) by about 30%. 

What could be the corresponding suppression in a strong-coupling scenario ? It is of 
course very difficult to answer this question given our impossibility to perform calculations 
in QCD at strong coupling. But if the experimental results at RHIC for the jet quenching 
parameter q |^9| , [30[ | — which, we recall, suggest an enhancement by roughly a factor of 5 with 
respect to the respective weak-coupling estimate — are indeed indicative of the strength of the 
quantum evolution in the QCD plasma, then one might expect a similarly strong reduction, 
by a factor of 5 or more, for the quark energy density in quenched QCD. That such an 
expectation is not totally unreasonable (within that strong-coupling scenario) can be also 
viewed via the following argument: 

Although there is no good reason to believe that the strong-coupling, large-n, estimates 
for the anomalous dimensions in J\f = 4 SYM theory, cf. Eq. (2.11), could be applied to the 
QCD problem at hand, let us nevertheless do so, by lack of a better argument. Previous 
studies in the literature, concerning the comparison between N = 4 SYM and thermal QCD 
in the temperature range of interest, suggest that a reasonable value for the QCD 't Hooft 
coupling to be used in this context is Aqcd — 5.5 [35, ^5|. (For instance, this is close to the 
naive estimate Aqcd = 3g 2 , with the 2-loop QCD running coupling g 2 (fl) evaluated at the 
scale Jl = 2ttT.) Via Eq. ( |2.11 ), this yields (for n = 2) an anomalous dimension |7//| ~ 1. 
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Assume now that there exists a window for strong-coupling dynamics, within which C^(/x 2 ) 
evolves according to Eq. ( |2,10| ), Then 

jfiV""^, (4.6) 



whereas the subsequent evolution from fi 2 to the harder scale Q 2 takes place at weak coupling, 
and hence it is much slower, cf. Eq. (|4.5|). Taking /xo = 3T C ~ 600 MeV once again, it is 
clear that a reduction by a factor of 5 or larger is achieved as soon as [i > 2/io ~ 1.2 GeV, 
that is, even if the strong-coupling dynamics holds only in a rather narrow window. The 
current lattice QCD results for the QCD pressure or energy density show a rather smooth 
behaviour for temperatures T > 3T C , with almost no variation from 3T C up to 6T C ; hence, if 
it so happens that QCD is (effectively) strongly-coupled at the scale 3T C , there is no reason 
why this should not remain true until the slightly harder scale of 6T C . 



To summarize, a lattice calculation for quenched QCD finding a result close to (4.3) 
would show that the "quasiparticles" of quenched QCD are close to being pointlike and that 
the theory is weakly coupled. On the other hand, a much smaller result, cf. ([4.2;), would be 
compelling evidence for an effectively strongly-coupled theory, with quasiparticles (if they 
exist) highly composite. 
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